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(D . Abstract 

We present parton-level analytical results for the next-to-leading order non- 
singlet virtual and real corrections to the Drell-Yan differential cross-section. 
The dependence of the differential cross section on the helicity of the ini- 
C^ I tial state partons is shown explicitly (the spins of the final state partons 

are summed). The calculation is implemented in dimensional regularization 
within the MS scheme and with the t'Hooft Veltman prescriptions for the 
n-dimensional 75. Both the polarized initial state and the unpolarized cross 
sections can be obtained from our result. Our unpolarized cross section agrees 
with the previous result of Ellis, Martinelli and Petronzio in the non-singlet 
sector. 
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1 Introduction 

Spin physics has emerged in the last few years as one of the most active fields of 
perturbative QCD. Given the planned experiments at RHIC and HERA, it is likely 
that such interest will continue to grow. Among the interesting processes which 
are planned to be investigated at hadron colliders such as RHIC are the Drell-Yan 
production of muon pairs and two-jet production in polarized proton-proton collisions 
[||. The difficulty of calculating spin dependence in next-to- leading order QCD is 
notorious and therefore only 0{as) results have been presented in the literature. 

The difficulty in the calculations comes from the presence of chiral fermions or 
gluons of fixed helicity in the initial state which renders the study of the spin de- 
pendencies extremely time consuming. The analysis of the infrared structure of the 
cross section in perturbation theory is also much more complex than the unpolarized 
case. If dimensional regularization is used to control both the infrared and ultraviolet 
singularities of the amplitudes, then the presence of chiral projectors for the fermions 
requires a suitable definition of 75 in n-dimension. 

In this work we present results for the O(a^) virtual corrections to the non- 
singlet sector of the classical Drell-Yan process. This is an important ingredient to 
the overall non-singlet differential cross section, da/dQ'^dtdu, and for the total cross 
section alQ"^), where Q^ is the invariant mass of the photon. The differential cross 
section is used to predict the transverse momentum distribution of the photon in the 
process p p —^ '~f* + X with polarized beams. 

It has been observed [^ that Drell-Yan production of charged lepton pairs is a case 
where transverse and longitudinal asymmetries are comparable. To leading order, 
the ratio between the transverse and the longitudinal asymmetries R = Att/^ll is 
estimated to be of order 1. This result is due to the fact that the gluons don't couple to 
transverse asymmetries and to order a^ only the quark-antiquark channel contribute. 
Furthermore, Drell-Yan is one of the possible avenues to obtain information about the 
twist-2 transversity distribution, which can be related to the ratio between transverse 
and longitudinal asymmetries. The inclusion of polarization in the initial states can 
be a useful tool to investigate such effects, and other effects, such as compositeness 
and supersymmetry (see and references therein). The presence of a massive photon 
in the final state makes an analytic calculation much more involved than, for instance, 
two other recent similar studies [Q. Although the non-singlet contribution to the 
process can only give information about the polarized quark distributions, it is the 
most lengthy Drell-Yan subprocess. We present here a brief account of the result for 
the order a^ corrections to the non singlet sector. 

2 General Structure 

The lowest order {non-singlet) contributions to the large transverse momentum pro- 
duction of virtual photons (see Fig. 1) with invariant mass q^ = Q^ arise from the 



two, q + q ^ 1* + 9, "Born" amplitudes shown in Fig. 3. We refer to these two 
diagrams as the direct and the exchange (or crossed) amphtudes, respectively. In 
Fig. 2 we have generically illustrated the expansion of the amplitude which appears 
in da/dQ'^ up to order a^. Notice that the quark form factor contributions repre- 
sented in Fig. 3 (and the related O(a^) corrections, not included in the picture) do 
not appear in the study of the cross section, da/d'^qtdy, for transverse momentum, 
qt, greater than zero. The two-loop 0{af) corrections to the quark form factor can 
be added in the study of the total cross section, da/dQ'^, by interfering the 2-loop 



on-shell quark form factor of ref. |]TI[ with the lowest order qq'j annihilation channel 
and by using helicity projectors for the initial quark states. 

We calculate the spin dependence of the cross section by using the helicity pro- 
jectors, P± = |(1 ±75), which project out the helicity states of an initial state quark 
and antiquark, respectively, as follows: 

m(P1,/Ii) = -(1 + /Ii75)m(p), V{p2,h2) = -V{p2){l + h2-f5), (l) 

where hi = ±1 = ±2Ai corresponds to quark helicity ±|, and /i2 = ±1 = ±2A2 
corresponds to antiquark helicity ±i. 

The squares of the direct amplitude Li and exchange amplitude L2 in Fig. 3 in 
N = 4 — 2e dimensions are given by 

M,,{h„ h2) = ej«^y ((1 - e)^ - hMl + ^r) , 

MUhu h2) = e)g'gl^- ((1 - ef - h^^ + ef) • (2) 

where a = g'^/'in is the fine structure constant, and cj the charge of the quark. 
The quantity Cp/Nc is the color factor, and a^ = gll^n is the QCD strong coupling 
constant. The interference term is more complicated and is given by. 



^^ ^' K. tu 



2Macih,h2) = 

- e)iQ^s - etu) -hih2{l + e){Q^s + etu) - 2hih2ttu . (3) 

The sum of the two Born amplitudes squared is 

|Mb(/ii,/i2)|' = \MB{h)\^ = Mrfrf(/i) +2M,,(/i) + M,,(/i) 

= ^)9'9l^^l^ [(1 - e) {2Q\s + (1 - e)(t^ + u^) - 2etu) 

-h{l + e) (2Q'^s + (1 + e){t^ + u^) + 2etu^ - Ahetu] , (4) 

and depends only on the product h = /ii/i2- 

To any order in perturbation theory we can write, 

\M{hi, /i2)r = \M{h)\^ = (1 - h)\M\^ + /i|M++|2 = |M|2 + hA\M\\ (5) 



where h = /11/12 = 4A1A2 and where 

1 

4 



hi,h2 hi,h2 



is the spin averaged (unpolarized) amphtude squared. Furthermore, 

|M__|2 = |M++p and |M_+|^ = |M+_|^ 

so that 

1 



The spin asymmetry, A|Mp, is defined according to 

1 



|Ma|' = - ( IM4 " 



IM, I 



M4 



|Mp 



(6) 

(7) 
(8) 

(9) 



The spin averaged (unpolarized) amphtude squared is determined from |M(/i)p 
by setting h = and A|Mp is the coefficient of h. For the Born term, this resuhs in 



\M\ 



2 22 2^F ^ 



1 - e) (2Q^s + (1 - e)(t2 + u^) - 2ttu 



and 



A|Mp 
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^/^ '^ N, tu 
Adding these two terms yields 



1 + e) [2Q'^s + (1 + e)(t^ + u^) + 2etM) + 4etM 



2 2 9 ^ F ^^ 



Q'^s + {t + uf 



(10) 



'111 



(12) 



which is proportional to e and vanishes in the limit e ^ 0. Since at the Born 
level there are no - singularities that might combine with this term to yield a finite 
contribution, in the limit e ^ 0, 



IM++P 



Oand A|M|2 = -|Mp 



(13) 



For the ? + ^ — ^ 7* + (? subprocess the condition that |M++p = means that the 
quark helicity is maintained (does not fiip) in the collision. The incoming quark line 
with helicity ±| turns around and becomes an incoming antiquark with helicity =f|. 
We refer to this as "helicity conservation". We see that at the Born level helicity is 
conserved in the limit e ^ 0. 

In A^ = 4 — 2e dimensions the differential cross section is related to the 2^2 
invariant amplitude according to 



d<y , ,x 1 /47rs\'' 1 ,,^/,m9 



dt 



I6ns \tu J r(l-e) 



(14) 



which can be written as 



or as 



da 1 \ da da^+ 

s-^(s, t, h) = {1- h) s—{s, t) + h s—^{s, t), 



where sda/dt is the unpolarized cross section and 



daiL _ 1 / da+_ 
^ dt ~ 2\ dt 



da. 



dt 



At the Born level we have, 



dt ■' s 

s^{s,t) = -e}K,^eAs{Q',u,t), 
dt •' s 



where 



TB{Q\u,t) 



2fl-e) 



1 - e) [2Q'^s + (1 - e){t^ + u^) - 2etu 
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tu - 

\t uj ut 



and 



where K2 is defined by 



AB(g^M,t) = — Q^s + {t + u) 

tu L 



Ko = na 



Cf 



2\ ^ / „/n2\ « 



47r/i^ \ I sQ 



N, r(l -e)\ Q^ 



tu 



(15) 
(16) 

(17) 



where we have rescaled, a^ — > ^^(/i^)^, so that it remains dimensionless in A^ 
dimensions. At this order of perturbation theory, we have 



s^(s, t, h) = e)K2— (1 - h)TB{Q'\ u, t) - heAniQ'', u, t) 



da , . , , o ^ , tt 

dt ' ■" s 

As we saw earlier, in the limit e ^ helicity is conserved so that 



da , , , da 



which implies that 
da. 



++, 



da 



LL 



da 



dt 



s,t) = and s^ — (s,t) = — s-— (s, t). 



dt 



dt 



(19) 
(20) 

(21) 
4-2e 

(22) 
(23) 
(24) 



To connect with the notation of ref . |^ , we note that 

K2TB{Q\u,t)=KTo{Q\u,t), 
where K and Tq are define in ref. [§] as 

Cf (1-e) [Atch^Y fsQ'^^ 



K = 2na 



N, r(i -t)\Q^ 



tu 



and 



^o(Q^n,t) 



t uj ut 



(25) 



(26) 



(27) 



3 Virtual diagram contribution 

The hst of the diagrams with virtual corrections contributing to the non-singlet sector 
of Drell-Yan is given in Fig. 4. We have omitted all the self-energy insertions of quarks 
and gluons and the ghost contributions. 

Our calculations are performed in the MS scheme using dimensional regulariza- 
tion to regulate both the ultraviolet and the infrared singularities. The isolation of 
the (unrenormalized) scalar invariant amplitudes at 1 loop level from the tensor inte- 
grals are done using the Passarino-Veltman procedure as implemented in FeynCalc 

i- 

We remove the ultraviolet singularities in the relevant subdiagrams by off shell 
regularization. Then by sending on shell the initial state quarks and the final state 
gluon, we encounter singularities in the form of double poles and single poles in e = 
2 — n/2. The one-loop order result for the cross section for the process q + (j —^ 'y* +g 
is given by 



da 



virtual 



dtdu 



[s,t,u, h) 



e}K2^5{s + t + u-Q')U{l- h)TB - heAe) 



1- 



a, r(l-e) (A-K^' 



27rr(l-2e) l Q^ 



.(H^^.i(3C.-.c.,n^,n,,,.,,£|_V 



a,9 



+ 2^(1-^) 



•.^(4C. + iVc) '^'':/^"' -2(2C, 



-2Cf 
-2 \Li2 



%{2Q^s + e + u^ 
tu 



3tu 
Q^s{t + u) 



N, 



c) 



Q^e + u" 
tu{t + u) 



t^ + u^ 



'rh^VH-'- 



tu{Q^ - u) (g2 _ t) (Q2 _ u) (g2 _ t) 



s^ + (s + u)'^^ 



^2lL,,|^^l+;i„Ml-^||liV,Hl±i^ + 2c/ + (^' + ')" 




t tu 

t + u\ 2Q^s + u^ + t^ + 2s'^ 



tu 



+2{2Cf - Nc) (lz, i^J-±^^ 

where h = hih2- The quantities Tb and Ab appear in at the Born level ( P^ and 
are given by equations (|19D and (|20|) , respectively. Notice that the double poles and 




single pole structure (in the /i — *> limit) automatically reproduces the singularities 
of the virtual contributions of ref. 0. In particular, the double poles in e correctly 
multiply the two-to-two Born contribution Tb, which is the Born level unpolarized 



result. Furthermore, it is possible to show |12] that the structure of the virtual 
contributions presented here correctly factorize, after adding the real contributions. 
The presence of the HsAb term in (^) implies that in the t'Hooft-Veltman |^, |T0| 
scheme the virtual corrections by themselves do not conserve helicity {i.e. they do 
not satisfy the equation (p3D). In particular, the finite part of daj^j^/dt in the limit 
e ^ for the virtual corrections is given by 

^^virtual 

o-^ / s 11 sO^ 1 \ 

-e]K2^^AB{s,t) h,CF - '^CF\n— + -Nc + Nc\n^- -Np] , (29) 

where s + t + u = Q^ . In some regularization schemes, such as those enforcing 
an anticommuting 75 in n dimensions, helicity would be conserved and equation 
(|13|) maintained throughout. In the t'Hooft-Veltman [^, ^ scheme helicity is not 
manifestly conserved for partial contributions which are not infrared safe. 



4 Real emissions 

Moving to the contributions from the real emissions in the non singlet sector in the 
same scheme, it is possible to show that they factorize according to the form 

[s,t,u, h) = 



dtdu 



a 



ejK2— I ((1 - h)TB - heAs) 6{s + t + u-Q' 



s 



as r(l — e) /47r/x 




+ (1 - h){dai + dan + dam) - ha^"^ - (1 + h)daiv, (30) 

where dai is the remaining finite part of the real qq annihilation diagrams denoted by 
dai in ref. [0], and dan, dam, and dajy correspond precisely to the finite unpolarized 
cross sections defined in [^. The term a^""^ is a regularization dependent term that 
arises from our use of the t'Hooft Veltman scheme. A more detailed study will be 



presented elsewhere |T^, however, the spin dependence of each term is shown. It 



is important to observe that the real term in eq. (pOf) contains the same factor of 
(1 — h)TQ — heAs as the virtual term in eq. (pSf). Except for the a'*"* term helicity 
is conserved along the quark (or antiquark) lines and the non-singlet order a^ spin 
dependent differential cross section can be written as follows: 

s-—-{s,t,h) = (1 - h){daj + dajj + dam) - (1 + h)dajv, (31) 

dtdu 

where daj, dau, dam, and dajy are the unpolarized cross sections given in ref. [0]. 
The regularization dependent term a^"'^ can be absorbed into the the spin dependent 
structure functions leaving a result that is exactly what is expected from helicity con- 
servation and a complete anticommuting scheme for 75. However, if this is done then 
one cannot use the MS scheme to calculate the Q^ evolution of the spin dependent 
structure functions (see ref. fT^). 



5 Conclusions 

We have presented the parton-level analytical results for the next-to-leading order 
non-singlet virtual plus real corrections to the Drell-Yan differential cross-section. 
Our unpolarized cross section, obtained from the h —>■ limit of (|31| ) agrees with the 
previous result of Ellis, Martinelli and Petronzio [Q] in the non-singlet sector, which 
we have recalculated. The dependence of the differential cross section on the helicity 
of the initial state partons is shown explicitly (the spins of the final state partons 
are summed). Although the calculation is very involved due to the presence of chiral 



projectors in the initial state, the result is quite compact and has been presented in 
a form from which cancelation of the infrared, coUinear and infrared plus collinear 
singularities is evident. Both the polarized initial state and the unpolarized cross 
sections can be obtained from our result. A more detailed discussion of our results 
and of the methods developed by us in the analysis of the virtual corrections will be 



presented elsewhere [12 



A study of the transversity distribution to O(a^) is also forthcoming [|T3 
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Figure 1: The Drell-Yan process 
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Figure 2: Radiative corrections of order a^ 
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Figure 3: The Born diagrams 
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Figure 4: One-Loop contributions 
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Figure 5: q + q-^'y* + G + G 
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Figure 6: g + g^7* + g + g 
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Figure 7: g + g-^7* + g + g 
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